Second order phase transitions induced by disorder in frustrated magnets 
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We study the critical properties of three dimensional frustrated magnets, diluted with non- 
magnetic impurities. We show that these systems exhibit a second order phase transition, cor- 
responding to a new universality class. In the pure case, the phase transition is expected to be 
weakly of first order. We therefore argue that these frustrated systems can be used to study experi- 
mentally the rounding effect of disorder on discontinuous phase transitions. We give first estimates 
of the critical exponents associated with this universality class, by using the method of the effective 
average action. 

PACS numbers: ll.10.Hi, 75.40.Cx 
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The influence of disorder on the properties of a system 
around a phase transition has triggered a lot of stud- 
ies in the last twenty years. In the case of second order 
phase transitions, we have a fairly good understanding of 
the situation thanks to the Harris criterion |[], 0], which 
states that for a broad class of systems with quenched 
disorder, the critical exponent v (describing the singu- 
larity of the correlation length at the phase transition) 
must satisfy the inequality: v > 2/d. This has a dra- 
matic consequence for these pure materials for which the 
preceding inequality is not satisfied, since a small amount 
of disorder changes the critical exponents, such that the 
preceding inequality is fulfilled. The pure system is there- 
fore unstable with respect to the inclusion of impurities, 
and one can expect a new universality class associated 
with the disordered system. 

Although much less studied, the influence of disorder 
on materials undergoing first order phase transitions is 
also of great interest. In such a situation, the general 
trend is that quenched disorder smoothens the phase 
transition. In particular, Aizenman and Wehr |3j have 
proven that when disorder is coupled to a system of small 
enough dimensionality p9| , the density conjugated to the 
disorder is a continuous function of the external param- 
eters (temperature, magnetic field, etc.). For a system 
with random bonds (i.e. a random distribution of the in- 
teraction between nearest neighbors degrees of freedom) , 
this statement implies that the energy density, which is 
then the density conjugated to the disorder, is a contin- 
uous function of the temperature. As a consequence, the 
latent heat vanishes: the transition is indeed smoothened 
by the disorder. 

This rounding effect of the disorder has been also stud- 
ied in specific models, mainly the random bond g-states 
Potts model. In two dimensions, Monte-Carlo simula- 
tions, as well as theoretical work show that the first order 
phase transition observed in the pure system for q > 5 
becomes of second order when quenched disorder is intro- 
duced (see 0] and references therein). Monte-Carlo sim- 
ulations on the three dimensional 3-states Potts model 
H display the same feature. However a shortcoming of 
the random bond Potts model for studying the round- 



ing effect of the disorder, is that it seems difficult to find 
experimental realizations in the relevant cases. Indeed, 
for three-dimensional realizations of the Potts model, the 
disorder usually induces a random field term, and thus 
the relevant model is not the random bond Potts model 
Q. On the other hand, in d = 2, it is required to find an 
experimental realization of the Potts model with at least 
five states, which seems quite difficult. 

While studying the rounding effect of the disorder, it 
would be valuable to find an experimental realization of 
this situation. With this goal in mind, we consider in 
this letter the influence of non-magnetic impurities in 
XY and Heisenberg stacked triangular antiferromagnets 
(STA), such as CsMnBr 3 , CsNiCl 3 , and rare earth he- 
limagnets (Ho, Tb, Dy). The spin degrees of freedom 
of these magnetic materials are submitted to competing 
interactions, which make the ground state more degener- 
ate than in a non-frustrated system. For the pure case, 
it was soon proposed that a new universality class should 
be associated with the phase transition observed in these 
materials |?J , and from then on, there has been numer- 
ous works aimed at characterizing it. We have now strong 
evidences that this phase transition is actually weakly of 
first order, and not continuous as was initially proposed. 
Let us indicate the most striking evidences for this behav- 
ior in the XY case. First, most experimental and numer- 
ical studies indicate power-law behaviors of the physical 
observables around the transition temperature, charac- 
terized by a set of critical exponents. A detailed analysis 
shows that these critical exponents are most probably not 
compatible with those expected in a second order phase 
transition || . It is more likely that the phase transition is 
actually of first order, with such a large correlation length 
that the discontinuity at the phase transition cannot be 
observed in experiments. This scenario is confirmed by 
a recent Monte-Carlo simulation || in which a first or- 
der phase transition was directly observed, but only for 
very large lattice sizes. Finally, a renormalization-group 
calculation in the framework of the effective average ac- 
tion shows a lack of a stable fixed point jlO| for both XY 
and Heisenberg cases, which in turn indicates that the 
phase transition is of first order. Note however that six- 
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loop resummed series indicate a true second order phase 
transition 

At the light of the introductory discussion, we propose 
to study the influence of quenched non-magnetic impuri- 
ties in the critical behavior of STA and helimagnets. In 
particular, we want to show that in presence of impuri- 
ties, the phase transition becomes continuous. We also 
want to characterize the associated universality class by 
giving a first estimate of the critical exponents. But be- 
fore embarking on this discussion, let us stress on two 
evidences which indicate the importance of disorder in 
these systems. The first observation is that, in the pure 
case, the effective exponent v (0.54 < v < 0.57 fL2| ) is 
found to be much smaller than the lower bound 2 /3 given 
by the Harris criterion [y < 2/d for disordered systems), 
so that we can expect the disorder to be strongly rel- 
evant. This conclusion is however an indication rather 
than a proof since Harris criterion only applies to second 
order phase transitions, which is not the case here. Sec- 
ond, the instability of the critical properties of the system 
when impurities are added has been tested experimen- 
tally on STA with diluted non-magnetic impurities [ jl3| , 
which shows, using scaling relation, that the exponent v 
increases strongly when impurities are introduced. 

We now derive the action relevant for the materials 
under study. Let us consider the general situation of 
a diluted system of spins on a lattice, submitted to a 
two-body interaction, in vanishing magnetic field. The 
corresponding Hamiltonian reads: 
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Here, Jjj represents the (translation invariant) interac- 
tion between spins at lattice sites i and j, and the spins 
Si are n-component vectors normed to unity. In physical 
realizations of STA and helimagnets, the number of spin 
components is taken to be two or three, depending on 
the anisotropies of the material. The dilution appears 
through the random variables q , which equal if a non- 
magnetic impurity lies on site i, and 1 otherwise. The 
Ci are chosen to be independent random variables with 
no dynamics, so as to describe a situation where the im- 
purities are quenched. The associated partition function 
reads: 



Z = / T>Sexp (-H) . 



(2) 



To deal with the unit norm constraint of Si, we perform 
the standard Hubbard-Stratonovitch transform (we fol- 
low the discussion of iQ) and introduce the auxiliary 
field if), so that the partition function reads: 

Z = J VSVi> exp -i J rj$Si + Yl 

(3) 



The variables Si at different lattice sites are now decou- 
pled, and the integral on S can be performed. One gets, 
up to an irrelevant multiplicative constant: 

Z = J IVexp f-i J2 J ^^3 - £ eiGM j (4) 



where [EOl 



^)=-in( /d "f£:i ) T^ ) ! (5) 



For deriving Eq.(Q) we used the fact that the q take 
values in {0, 1}. In terms of the local magnetization tpi = 
J2j Jrl^j @' ^ e partition function reads: 



(6) 

The important point here is that G n (x) is a scalar un- 
der rotations (see Eq. (|g)) and therefore depends only 
on x 2 . We conclude that the disorder ej couples to the 
microscopic field ip through terms with even powers of (p 
such as, for instance, the random mass term e^f , but not 
through terms linear in the field, which would correspond 
to the more complex situation of a random field. 

We now specialize to the case of helimagnets and STA. 
In these materials, the spin degrees of freedom are sub- 
mitted to competing interactions, which induce frustra- 
tion. As a consequence, there are two critical modes 
L P±Q (*' e ' m °des minimizing the two-point interaction 
J(q) in Fourier space), with Q some characteristic mo- 
mentum H instead of only one critical mode (/%=o m 
the ferromagnetic case. In the long distance limit, the 
physics is governed by those modes whose momenta are 
very close to ±Q, while the other modes only influence 
non-universal quantities, and can be discarded. We are 
therefore lent to introduce two types of fields: 



<p±(o) = v^Q + q)- 



(7) 



(8) 



It is convenient to perform the change of variables: 

<pi(q) = <p + (q) + i<p_(q) 

<P2(q) = <p+(q) - i<P-(q) 
If we now expand the interaction potential G n to fourth 
order in the fields, the partition function can be recast 
into the form: 



Z = J T>0iT><p2 exp(— S) 



(9) 



with: 



S = 



flips) 2 



-(x) 



92{x) 



(10) 
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We have merged here the two n-components vectors (p\ 
and Cp 2 into a 2 x n matrix <fi = (</5i, (p 2 ), and introduced 
p = Tr (*</>. 4>), t = Tr This action is very 

similar to the one obtained in the pure case (7) , except 
for the x-dependence of the field normalization Z and of 
the coupling constants r, gi and g 2 , appearing through 
the presence of impurities (noted e$ in Eq.(|l|)). In the 
following, we consider a direct generalization of this situ- 
ation where </> is a m x n matrix. For to — 2, we retrieve 
the case of physical interest discussed previously. 



exact, concise equation Ifi 



Once the action (10) is determined, one would have to 
compute the quantity of interest for a given realization of 
the disorder (for a particular choice of the {ei}). This is 
however a formidable task. A simpler strategy consists in 
studying self-averaging quantities averaged over the dis- 
order. By doing so, one recovers translational invariance. 
We use here the replica trick for computing the average 
of the free energy over the disorder. We write: 



F = In Z = lim 

0^0 



Z° - 1 



(11) 



where the line over a quantity represents the average with 
respect to the disorder. The o th power of the partition 
function can be written by introducing o copies of the 
field 4>. The average over disorder can now be computed 
(at least for a Gaussian noise), and if we only keep the 
marginal and relevant terms in the sense of renormaliza- 
tion group, we find the following "0 4 -like" action: 



S ^ k ) = |d rf J]T[§Tr(9 
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where the index I labels the o copies of the field. The 
critical properties of the materials under study are ob- 
tained by using the previous action, and taking the limit 
o — > at the end of the calculation. 

Various methods can now be used for studying the crit- 
ical properties of the system (expansion in coupling con- 
stants, e-expansion, etc.). We propose here an approach 
based on the concept of effective average action. This 
is a well documented topic (see jl5| for a review, and 
JTo| , 16 for applications in closely related problems), so 
we only present the main ideas here. This approach is 
based on the study of the effective average action Tk, 
which interpolates smoothly between the microscopic ac- 
tion S (when k — > 00) and the Gibbs free energy T (when 
k — ► 0). Tk therefore connects smoothly the macroscopic 
and microscopic descriptions of the system. For interme- 
diates values of A:, Tj~ is obtained by integrating over the 
high momentum modes only, while long distance modes 
are effectively left unintegrated. A major interest of Tk 
is that its evolution with the scale k is governed by an 



dTk 1 f d t R k (q) 
8t 2j q m(q)+R k (qY 



(13) 



with k oc exp(i). The function Rk{q) describes how the 

(2) 

low and high momentum modes are separated, and T k 
is the second functional derivative of the effective average 
action. Although this equation cannot be solved exactly, 
it can be used in practice to determine nonperturbative 
flow equations. The idea consists in considering a trun- 
cation for r^, typically characterized by a finite number 
of coupling constants. By introducing this truncation 
in equation (13), one deduces flow equations for the dif- 
ferent coupling constants. These are non-polynomial in 
coupling constants, as can be seen from the nonlinear 
structure of the flow equation (13). 

We now describe the type of truncations considered 
here: 



r * = / f E Tr ( d + y E to - ^ 



+ f!>-^) + f (e 



(14) 



(pl - «) 



This truncation is to be seen as the first terms of an 
expansion of the effective average action around a con- 
figuration M in, defined by {<t>uxni,j)l — V^fiij (here, i 
and j, running respectively from 1 to n and from 1 to to, 
stand for the matrix indices of the field, while I stands 
for the replica index, running from 1 too). We have also 
considered truncations with <j) e terms added. The action 
depends then on six more coupling constants. The flow 
equations have been derived using the method presented 
above. However those are to lengthy to be displayed here. 
We only quote the flow equation in the limit of small cou- 
pling constants {ui, 1*2,1*3}, which identify with the one 
loop expansion in coupling constant [pl| : 



u\ = —sui 



+ 2v 4 ((8 + mn)u\ + V2u\ + (15) 
+4ui((l + m + n)u 2 + 3u 3 )J 
«2 = -eu 2 + 4w 4 m 2 ( 6(mi + u 3 ) + (4 + to + n)u 2 ^j (16) 
U3 = — eu3+2v4U3 ^2(2 + mn)u\ + (17) 
+ 4(1 + to + n)u2 + (8 + mno)u^j 

where e = 4 — d, V4 = 1/327T 2 , and the dot indicates a 
derivative with respect to t. For generic values of m, n 
and o, the flow equations admit eight fixed points. Out 
of them, six are already well known. There are four fixed 
points with u 2 = 0, corresponding to the situation of 
a diluted system of mn components spins |lS| ], and two 
extra fixed points with 113 = 0, corresponding to the pure 
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frustrated magnets |6|. In addition, there are two new 
fixed points. We finally indicate that for certain values of 
(to, n, o), there is a degeneracy of the one loop (3 functions 
(for example, for to = 2, o = 0, n ~ 25.3). For such 
values of (to, n, o), we expect a -y/i-epsilon fixed point to 
appear at two-loop order, as in the random bond Ising 
model Q. We postpone the discussion of this point, as 
well as a more detailed description of the fixed points to 
a forthcoming publication. 

In three dimensions, we have looked for roots of the j3 
functions. In the case of physical interest m = 2, the dis- 
order is found to be irrelevant for n > 5.5. On the other 
hand, for n < 5.5, the disordered fixed point is stable 
and governs the phase transition. We give on Figjl] the 
values of the associated critical exponents. We first re- 



particular, a new universality class was shown to appear, 
with associated critical exponents differing strongly from 
the effective exponents measured in the pure case. We 
also gave some indications for a very appealing behav- 
ior of the critical exponent v, and a large value of r\. A 
verification of these effects by independant approaches, 
in particular perturbative expansions and Monte-Carlo 
simulations, would be very valuable. At the light of the 
results presented here, it would be particularly interest- 
ing to reconsider the experimental studies on diluted frus- 
trated magnets. 

We thank B. Delamotte, D. Mouhanna and M. Picco 
for useful discussions. This work was supported by Marie 
Curie fellowship HPMF-CT-2001-01343. 
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FIG. 1: Critical exponents of the diluted STA, as a function of 
the number of spin components. The dotted line corresponds 
to the lower bound given by Harris criterion. 

mark that our values of v are all smaller than the lower 
bound 2/3 predicted by the Harris criterion. However, 
for 3<n<5.5,z^is off by only few percents, which is 
typically the error expected with the truncation consid- 
ered here. The much larger discrepancy for n — 2 can be 
attributed to the unusually large value of the anomalous 
dimension to In this case, richer truncations are to be 
considered in order to get more accurate results (work in 
progress) . Although the accuracy of our results is ques- 
tionable, in particular for n < 3, they indicate that v is 
varying very slowly with n (in the pure vectorial model, 
the variation of v on the same range of n is almost ten 
times bigger), and that v is close to the lower bound 2/3 
imposed by Harris criterion. These features are similar 
to what is observed in d = 2 gr-state Potts model with 
random bond using Monte-Carlo simulations, where v is 
found to be very close to the lower bound 1 = 2/d ob- 
tained from the Harris criterion, for all q > 2. 

To conclude, we have highlighted some interesting 
properties associated with the inclusion of non-magnetic 
impurities in three-dimensional frustrated magnets. In 
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